Under consideration for publication in J. Fluid Mech. 



1 



Inertial effects in the SafFman- Taylor 

instability 

By CHRISTOPHE CHEVALIERi^, 
MARTINE BEN AMAR^, DANIEL BONN^'t 
AND ANKE LINDNERi 

""^Laboratoire de Physique et Mecanique des Milieux Heterogenes, 
UMR 7636 CNRS, Ecole Superieure de Physique et de Chimie Industrielles, 
10 rue Vauquehn, 75231 Paris cedex 05, France 

^Laboratoire de Physique Statistique, UMR 8550 CNRS, Ecole Normale Superieure, 

24 rue Lhomond, 75231 Paris cedex 05, France 
^Ecole Nationale des Fonts et Chaussees, 6-8 avenue Blaise Pascal, Cite Descartes, 
Champs sur Marne, 77455 Marne-la-Vallee cedex 2, France 

*Van der Waals-Zeeman Institute, University of Amsterdam, 
Valckenierstraat 65, 1018 XE Amsterdam, the Netherlands 

(Received 2 February 2008) 

For the SafFman- Taylor instability, the inertia of the fluid may become important for 
large Reynolds numbers Re. We investigate the effects of inertia on the width of the 
viscous fingers experimentally. We find that, due to inertia, the finger width can increase 
with increasing speed, contrary to what happens at small Re. We find that inertial effects 
need to be considered above a critical Weber number We. In this case it can be shown 
that the finger width is governed by a balance between viscous forces and inertia. This 
allows us to define a modified control parameter 1/B\ which takes the corrections due 
to inertia into account; rescaling the experimental data with 1/5', they all collapse onto 
the universal curve for the classical Saffman- Taylor instability. Subsequently, we try and 
rationalize our observations. Numerical simulations taking into account a modification of 
Darcy law to include inertia, are found to only qualitatively reproduce the experimental 
findings, pointing to the importance of three-dimensional effects. 



1. Introduction 

Viscous fingering has received much attention as an archetype of pattern-formation 
problems and as a limiting factor in the recovery of crude oil fsee ISaffman fc Tavloil 
ll958HBensimon et oiJlT98ft lHomsvlll987t ICoud'alllOQll) . Viscous fingers form when in a 
thin linear channel or Hele-Shaw cell, a fluid pushes a more viscous fluid. The interface 
between the fluids develops an instability leading to the formation of finger-like patterns. 
The viscous fingering instability has been studied intensively over the past few decades 
both theoretically and experimentally. 

For the classical Saffman- Taylor instability the width of the finger is governed by the 
competition between viscous and capillary forces: viscous forces tend to narrow the finger 
whereas capillary forces tend to widen it. When air pushes a viscous fiuid, as is usually 
the case, the relative finger width is thus determined by the capillary number Ca = r/U/j, 
(with 77 the viscosity, U the velocity and 7 the surface tension) the ratio between viscous 
and capillary forces. In the vast majority of cases that have been studied so far, inertial 
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forces are negligible. The importance of inertia is given by the relative importance of the 
inertial and viscous forces, quantified by the Reynolds number Re = pUb/r], with p the 
fluid density, and b the plate spacing of the Hele-Shaw cell in which the experiments are 
conducted. In most studies of the instability, b is small, and the fluids considered both 
in applications as well as in experimental studies are typically high-viscosity oils. This 
automatically leads to small Reynolds numbers {Re << 1), so that inertial effects may 
be neglected. 

More recently viscous fing ering has been studied in non-Newtonian fluids using for ex- 
ample polymer solutions (seelSmith et aZ.ll993lBonn gJlQQStlLindner. Bonn &: Meunieil 
l2nnnt IVlad fc Mah'CTll200ft iKawaguchi. Hibino fc Katoll200lt iLindner et QiJl2002h . For 

the dilute polymer solutions used in a number of these studies, the shear viscosity of 
the water-based solutions is typically close to the water viscosity, and consequently the 
Reynolds number may -and does- become larger than unity. This means that inertia may 
become important, and needs to be disentangled from the observed non-Newtonian flow 
effects. Also recently, correctio ns to the so called Darcy ' s law have been developed in- 
corporating inertial effects fsee lOondret fc R,abaudlll99^lR,uver-Qui]lEoOll) . Darcy's law 
relates the pressure gradient to the fluid velocity and is one of the fundamental equa- 
tions of the Saffman- Taylor instability; if inertial corrections can simply be included in 
a modified Darcy's law, this would greatly facilitate the understanding of the effect of 
inertia on the instability. These recent developments suggest that a better understanding 
of the Newtonian fingering instability for high Reynolds numbers is both necessary and 
feasible. 

In this paper we explore the Saffman- Taylor instability for Newtonian fiuids for Reynolds 
numbers up to Re = 100. To do so, we use low- viscosity silicone oils, pushed by air. 
The paper is organized as follows. In Sec. 2 we will recall the basic equations for the 
Saffman- Taylor instability and introduce the corrections due to inertia. Sec. 3 describes 
the set-up and experimental methods. In Sec. 4 the experimental results concerning the 
finger widths as well as the validity of Darcy's law are presented and discussed. In Sec. 5 
wc will introduce some theoretical elements as well as numerical simulation and their 
comparison to the experimental results. Sec. 6 gives a summary of the obtained results. 



2. Theory and equations 

2.1. Presentation and review of classical Saffman- Taylor instability 

We study the Saffman- Taylor instability in a thin linear channel or Hele-Shaw cell (see 
figure^. The width of the cell W is chosen to be large compared to the channel thickness 
b and we thus work with high aspect ratios W/b. The cell is filled with a viscous fluid 
which is subsequently pushed by air. The viscosity and the density of air will be neglected 
throughout the paper. 

When air pushes the viscous fluid due to an imposed pressure gradient VP, an ini- 
tially flat interface between the two fluids destabilizes. This destabilisation leads to the 
formation of a so called viscous finger; in steady state a stationary finger of width w 
propagating at a velocity U is found to occupy a fraction of the cell width: the relative 
finger width is defined as: A = w/W. 

For Newtonian fluids, the motion of a fluid in the Hele-Shaw cell is described by the 
two-dimensional velocity field u averaged through the thickness of the cell. It is given by 
Darcy's law, which relates the local pressure gradient to the velocity within the fiuid as: 

u^-^Vp. (2.1) 
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Figure 1. Schematic drawing of the experimental set-up. 



It follows immediately that, if the fluid is incompressible, the pressure field satisfies 
Laplace's equation: 

Ap = 0. (2.2) 

The pressure field is calculated within the driven fluid with in addition a pressure jump 
over the interface due to the surface tension: 



5p = 7/i?, 



(2.3) 



with R the radius of curvature of the interface, again employing a two- dimensional ap- 
proxim atio n, as was justified in the l imit of small capillary numbers bv IPark fc Homsvl 
lll984l) . and lReinelt fc SaffmanI (119851) . 

The other boundary conditions are the continuity condition, which implies that the 
normal velocity at both sides of the interface is equal and a far-field value for the pres- 
sure. Supplemented with these boundary conditions, H2.1|l . (|2.2(l and H2.3|l constitute the 
complete set of equations one can solve in order to obtain the complete finger shape for 
a given pressure gradient; and thus also its width. 

For characterizing the instability quantitatively most studies have focused on the width 
of the finger w relative to the channel width W: A = w/W as a function of their velocity. 
It follows from the above that their width is determined by the capillary number; one 
thus anticipates that the relative width of the viscous fingers decreases with increasing 
finger velocity. This is indeed what is generally observed experimentally; in addition, 
for very large values of Ca, A does not go to zero but reaches a limiting value of about 
half the channel width. It also follows from the boundary conditions and (|2.1() to (|2.3() . 
that the control parameter for the fingering problem is = 12{W/b)^ Ca with W/b 
the aspect ratio of the Hclc-Shaw cell. When scaled on 1/B, measurements of A for 
different systems all fall on the same universal curve. In the ideal, Newtonian, two- 
dimensional situation 1/B is consequently the only parameter th at determines the fin- 
ger width fsee .Saffman fc Tavlor 1958: ..McLean fc Saffman,. 1 98 L: iCombescot et a?j|l986l: 
iHong fc Langeilll98(iHShraimanlll98(f . 

2.2. Corrections due to inertia 

When inertial forces have to be taken into account, both the Reynolds number Re = 
pUb/rj or the Weber number We = pU'^b/^ (the ratio of inertial forces to capillarity 
forces) may become important. We will now discuss how corrections due to inertia can 
be included in the basic equations. 

Modifications of Darcy's law have first been proposed by iGondret fc RabaudI l)l997() 
for parallel flow in a Hele-Shaw cell: they establish corrections by averaging inertia in 
the third dimension, i.e., they average over the direction of the plate spacin g 6, allowing 
them t o derive a new nonlinear two-dimensional equation for the velocity field. iR.uver-Quill 
suggests an improved correction starting from the three-dimensional Navier-Stokes 
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equation. Incrtial corrections arc introduced in a perturbative fashion; using in addition 
a polynomial approximation to the velocity field, Ruyer-Quil proposes the modified two- 
dimensional Darcy's law of the form: 

with a = 6/5 and /? = 54/35. IPlouraboue fc HinchI ll2002l) also calculated inertial correc- 
tions to Darcy's law and arrive at a similar type of equation, but with slightly different 
coefficients. This equation leads to a better agreement between the linear stability analy- 
sis and the experimental data of Gondret & Rabaud for the Kelvin-Helmholtz instability 
up to not too large Reynolds numbers, a and (3 may vary depending on the way the 
averaging in the third dimension is done, but are always of order of 1. 

When scaling length on VF, time on W/U and pressure on 12rjUW/lP' one finds the 
following dimensionless equation: 

i.e*(.^+/5«*V*n*)^-VV-«*, where R^* = y,^'-^ = ^R^- 

(2.5) 

Re* is a modified Reynolds number, in the same way as the classical control parameter 
of the Saffman- Taylor instability 1/5 is a modified capillary number. 

We can also introduce another number describing the importance between inertia and 
capillarity in the geometry of the Hele-Shaw, a modified Weber number: 

We* ^ f^^-— = -^We. (2.6) 

7 

One important remark is that if one considers stationary and spatially uniform flow in 
our Hele-Shaw cell, it follows from (|2.4() that there are no corrections due to inertia, since 
du/dt and u^u are both zero. This will be the case in our fingering experiments far 
away from the moving interface and leads to the classical Darcy's law; we thus anticipate 
that it might remain valid even for relatively high Re. 



3. Experimental 

We use a linear Hele-Shaw cell consisting of two glass plates separated by a thin Mylar 
spacer. The plates are horizontal and clamped together in order to obtain a regular 
thickness b of the channel. The thickness of the glass plates is chosen to be 2 cm in order 
to avoid any bending of the plates. The aspect ratio of the channel can be varied; we 
worked with different plate spacings b and widths of the cell W , the length of the channel 
being always 1 m. The cell is filled with silicone oil and compressed air is used as the less 
viscous driving fluid. 

The silicone oils used were RhodorsU 47V05, 47V10, 47V20 and 47V100 from Rho- 
dia Silicones. Rheological measurements on a Reologica Stress- Tech rheometer confirmed 
the values of the viscosities of 5, 10, 20 and 100 mPa.s respectively, with no devia- 
tions larger than 4%. We also used 47V02, its viscosity was measured to be 2.8 mPa.s. 
The surface tension 7 and the density p of the silicone oils are 19.5 ±1 mNm^^ and 
0.95 ±0.03 10"'^ kgm~'^ as given by Rhodia Silicones. 

The fingers were driven by applying a constant pressure drop Ap — Pi — Po between 
the inlet and the outlet of the cell. Depending on the order of magnitude of the applied 
pressure drop two methods were used. For Ap larger than 3000 Pa we used compressed 
air and a pressure transducer at the entrance of the cell to fix pi at the inlet of the cell. 
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Thickness b Width W 



Geometry 1 
Geometry 2 
Geometry 3 
Geometry 4 



0.25 mm 
0.75 mm 
0.75 mm 
1.43 mm 



40 mm 
80 mm 
40 mm 
40 mm 



Table 1. Different cell geometries used in our experiments. 



In this case, the outlet was maintained at atmospheric pressure Po = Patm due to an oil 
reservoir coupled to the cell. For a Ap smaller than 3000 Pa, we obtained the pressure 
drop by lowering the oil reservoir at the outlet of the cell by a given amount, determining 
in this way Po- In this case the inlet was maintained at atmospheric pressure pi = Patm- 

The fingers were captured by a CCD camera, coupled to a data acquisition card (Na- 
tional Instruments) and a computer. This allowed for measurements of the relative width 
A = w/W as a function of the velocity U of the finger tip. For each configuration (cell ge- 
ometry and fluid viscosity) several experimental runs (between 10 and 20) were performed 
increasing the applied pressure drop and thus the finger velocity until destabilization of 
the finger occurred; all the finger widths reported here correspond to stable fingers. 

In order to access high Reynolds numbers we can not only vary the velocity of the 
finger and the viscosity of the fluid but also change the thickness of the channel. We have 
thus worked with different channel geometries that are summarized on tabled 

For the geometries used, the aspect ratio W/b varies from 28 (geometry 4) to 160 
(geometry 1). Even if an aspect ratio of 28 is rather small it is sufficient to consider the 
experiment as being quasi two-dimensional. It is observed that the results obtained for 
the high viscosity fluids (and thus a situation where inertial effects can be neglected) 
show only yerv Httle difference in their finger widths. The small difference is due to film 
effects fsee iTabeling fc Libchabeijll98d) . as will be discussed in more detail below. 

Experiments were performed in all geometries for the silicon oils 47V05, 47V10 and 
47V20. The Silicon oil 47V02 was tested in geometries 2 and 3, whereas the silicon oil 
47V100 was used in geometries 3 and 4. Finally, note that typical values of capillary 
number Ca in the experiments are between 0.01 (for V02) and 0.5 (for VlOO). 



4. Presentation of the results 

4.1. Darcy's law 

Assuming the flow far away from the finger to be uniform, we expect that the classical 
Darcy's law linking the gap averaged fluid velocity V to the imposed pressure gradient 
VP in our Hele-Shaw cell remains valid for all of our experiments: 

V = ~^VP. (4.1) 

Mass conservation allows to obtain the velocity V of the fluid far away from the 
interface from the measured finger velocity U simply by using V = XU, if one neglects 
the thin wetting film left on the glass plates behind the finger. The imposed pressure 
gradient is calculated by: VP = Ap/L where Ap is the measured applied pressure drop 
and L the distance between the finger tip and the exit of the cell. 

In our experiments we reach high finger velocities and thus high capillary numbers Ca. 
The infiuence of the thin wetting film left on the plates may therefore become important 
and can not be neglected any more in some of the experiments. It is taken into account 
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Figure 2. Velocity as a function of the applied pressure gradient for all viscosities (V02, V05, 
VIO, V20 and VlOO) and different cell geometries: •, 6=1.43 mm VK=4 cm; x, 6=0.75 mm 
W=A cm; □, 6=0.75 mm W=8 cm; +, 6=0.25 mm W=4 cm. 

using V = XU{1 — 2t/b), wh ere t is the thickness of the w ettinK film, which we estim ate 
using the empirical result of IXabeling fc Libchabeil ||1986^ and Ta bcling et all l|l987|) : 

t = Kb[l - cxp(-7iy/6)][l - exp(-/3Co^/^)], (4.2) 

with K « 0.119, 7 « 0.038 and f3 w 8.58. 

In this way, we can thus test the validity of Darcy's law. Figure |2| shows the velocity 
V represented as a function of {b/12r])VP for the different cell geometries and viscosities 
used. The dashed line represents the linear relation with slope unity expected from H4.1|l . 
We therefore conclude that the data are in excellent agreement with the classical Darcy's 
law. This result holds even for high velocities where a significant effect of the inertial 
forces is observed on the width of the fingers, as will be discussed below. We have thus 
shown that for the range of Re tested in this paper there is, as was anticipated above, 
no effect of inertia on Darcy's law when considering the uniform flow far away from the 
finger. Note that this does of course not mean that there are no corrections to the local 
Darcy's law near the finger tip. 

4.2. Finger width 

4.2.1. Relative finger width as a function of the classical control parameter 

Figures |2Ia) andl^b) represent the relative finger width as a function of the classical 
control parameter 1/B when varying the viscosity of the fluid for a given geometry - 
6=0.75 mm W=A cm, figurelHIa)- and when changing the geometry of the cell for a given 
fluid -silicon oil 47V05, figureEIb). These figures show that, for low 1/B one observes the 
classical decrease of the finger width with increasing 1/B. However at a given value of 
1/B which is different for different configurations, an increase of the relative finger width 
is observed. This surprising observation systematically appears at high Reynolds number, 
and we conclude that it must be related to inertial effects. Indeed, for a given geometry, 
only the fluid of highest viscosity gives results that agree with the classical Saffman- 
Taylor instability. In addition, deviations from the classical results arise at smaller 1/B 
for lower fluid viscosity. Finally, the data for a fixed viscosity but varying geometry (figure 
EJ)) show that the increase of the finger width occurs for lower 1/B for a thicker channel. 
All these observations agree with the suggestion that the increase in finger width with 
increasing velocity is due to inertial effects. 
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Figure 3. Results for the finger width A as a function of the classical control parameter 1/B: 
(a) for the geometry with 6=0.75 mm and W=4 cm and different fluids: •, V02; x, V05; □, VIO; 
+ V20; A, VIOO. (b) for the V05 fluid and different cell geometries: •, 6=i.43 mm W^4 cm; 
X, 6=0.75 mm W=4: cm; □, 6=0.75 mm W=8 cm; +, 6=0.25 mm W=4 cm. 

Comparing the data for a fixed gap thickness (5=0.75 mm) and two different channel 
widths (Vl^=4 and W=8 cm) on figure |3Jb), it follows that the crossover value of 1/B 
also depends on the channel width W: it is observed to be smaller for smaller channel 
width. This is still consistent with an increase of the Reynolds numbers {Re and Re*): 
at a given 1/B and fixed rj and b, a decrease of the channel width W leads to an increase 
of Re. This follows from the observation that W^U is fixed and consequently Re varies 
as 1/W (Re* as 1/W^). 

We conclude that due to inertia our experimental results deviate from the classical 
curves: we observe a regime of increasing finger width. This increase occurs at lower 1/B 
for lower viscosity, larger gap thickness or smaller gap width. It is also important to 
note that strong inertial effects are observed already at velocities below 0.08 ms~^ for 
all of our experiments and that even if strong changes in the behaviour are observed for 
the finger width, no deviations from the classical Darcy law is observed in this regime 
(figure 121). 

We will now study the deviation from the classical result as a function of the modi- 
fied Reynolds number Re*. For concreteness, in the following we will focus on the data 
obtained when varying the viscosity for a given geometry (5=0.75 mm, W=A cm). The 
results however are general and apply to the data from other experiments as well. 

4.2.2. Relative finger width as function of the modified Reynolds number Re* 

In figure ^a), we plot the relative finger width as a function of the modified Reynolds 
number Re*. The first observation is that the minimum of the curves that signals the 
deviation from the classical results, is not given by Re* . On the other hand for high Re* 
all the curves tend towards a single master curve: the behaviour of the finger width seems 
to be governed by Re* only. Data in other configurations (not shown here) confirm the 
existence of a universal A-i?e*curve. 

4.2.3. Relative finger width as a function of the modified Weber number We* 

So far we can distinguish between two limiting cases. For low velocities, the results 
for the relative finger width fall on the universal curve of the classical Saffman- Taylor 
instability: they rescale with 1/B. For high values velocity, a second universal curve exists 
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Figure 4. Results for the finger width A as a function of the modified Reynolds number Re* (a) 
and of the modified Weber number We* (b), for the geometry with 6=0.75 mm and W—4 cm 
and different fluids: •, V02; x, V05; □, VfO; + V20; A, VfOO. 



and the data rescale with Re*. It follows that the crossover between the two regimes is 
given by the modified Weber number, combining Re* and 1/B: 

pU'^W W 

We* = Re*. 1/B = ^ = — We. (4.3) 

7 b 

The experimental data support this conclusion. Figure 2fb) depicts the relative finger 
width as a function of the modified Weber number We* . All experimental curves have 
a minimum located at the same value of We* at around We^ « 15 separating the two 
limiting behaviours. 

Note that, although We* governs the crossover, we observe no regime where the fin- 
ger width is given by a competition between capillary forces and inertia. In fact, when 
considering the dependence of the different forces on the velocity one finds that capillary 
forces scale as t/", viscous forces as and inertial forces as t/^. Consequently, the dom- 
inating forces should at low velocity be capillary and viscous forces (control parameter 
1/B) and at high velocity, viscous forces versus inertia (control parameter Re*). This 
simple argument therefore explains that as a function of the velocity there is no regime 
where the finger width is given by We*. 

4.2. A. Extension to a new global master curve 

It follows that the parameter We* {— Re*. 1/B), that can be seen as the ratio between 
1/B and 1/Re*, gives the relative importance of the two parameters with a cross over 
given by the critical value We* « 15. 

We can thus attempt to define a modified control parameter taking this crossover into 
account: 

1/B' = 1/B ( —, V (4.4) 

' ' \1+We*/We*,j ^ ' 

It is easily seen that this parameter tends to l/S for low We* {We*< We*) and towards 
We*JRe* for large We* {We*>We*). 

Figure shows the experimental data already shown on figure |2Ia), however now A is 
plotted as a function of 1/5'. Surprisingly the experimental data scale on a single uni- 
versal curve when represented as a function of the modified control parameter. Moreover, 
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Figure 5. Results for the finger width A (same data as figure]^) as a function of the modified 
control parameter 1/B\ for the geometry with 6=0.75 mm and 14^=4 cm and different fluids: 
• , V02; X, V05; □, VfO; + V20; A, VfOO. [Insert: large scale] 
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Figure 6. Results for the finger width A as a function of the classical control parameter 1/ B (a) 
and of the modified control parameter (b) for aU viscosities (V02, V05, VIO, V20 and VlOO) 
and difi'erent cell geometries: •, 6=1.43 mm W=4: cm; x, 6=0.75 mm W=4: cm; □, 6=0.75 mm 
VK=8 cm; +, 6=0.25 mm W=4: cm. 



and perhaps even more surprisingly this curve is identical to the result of the classical 
Saffman- Taylor instability. 

The figures E{ a) andlHfb) summarize our results. On figure |^a), all geometries and 
fluids are depicted representing the relative finger width as a function of the classical 
control parameter 1/B. One should emphasize that this data are obtained by varying 
not only the fluid viscosity but also the cell geometry by changing both, the channel 
width and thickness. On figurcEIb), the same data are plotted as a function of 1/B', our 
modified control parameter. We observe that the entire data set collapses very reasonably 
onto the single master curve when represented as a function of 1 /B\ Once again this 
curve is identical to the classical result of ISaffman fc Tavloil l)l958j) . 

So far we did not discuss the infiuence of the aspect ratio on the relative finger width. 
Even if the infiuence of the latter is small it might explain the fact we do observe slight 
deviations from the four different cell geometries (see figure Eb). However when consid- 
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Figure 7. Snapshot of a finger: (a) without inertial effects {We*— 5 < Wei, silicon oil V02, 
velocity f/=5.2 cms"'^, 6=0.75 mm, W=4: cm); (b) with inertial effects ( VFe*=30 > Wei, silicon 
oil V02, velocity ?7=f2.5 cms "*^, 6=0.75 m m, W=4cm). x: Finger shape predicted without 
inertial effects by the theory of lPittj (If 98011 . 



ering one single channel geometry the data do collapse (see figure O. If one also notes 
that these differences can already be observed where inertia is negligible, the conclusion 
must be that the slight residual differences are due to film effects. 

The physical interpretation of our results is then the following. The modified control 
parameter 1/B' gives the cross-over between 1/B and Re*. For small We*, 1/B is the 
control parameter, and the main forces are surface tension and viscous forces, leading 
to a narrowing of the fingers width as viscous forces become more important for higher 
speeds. For higher velocities {We*>We*), the main acting forces become viscous forces 
and inertia. The competition between these forces results in a widening of the finger 
width with increasing velocity. The observation is therefore that inertia tends to widen 
the fingers; this seems logically intuitively, as the inertia will tend to slow down the finger 
at a given flow rate, leading consequently to wider fingers. As the effect of the inertial 
forces is similar to that of the capillary forces in, the sense that both tend to widen the 
finger, and the classical Saffman- Taylor finger selection appears to have remained intact, 
one may attempt to include the inertial forces in an effective surface tension. Indeed the 
modified control parameter 1/B' can be written as the classical control parameter 1/B 
by including an effective surface that is of the form: 

^eff =1(1+ We*/We*J. (4.5) 



5. Some theoretical elements 

In addition to the above, the experimental observations indicate that even if the finger 
width increases when increasing the velocity sufficiently the finger shape does not really 
change (see figure [TJ . All these observations suggest the possibility to introduce the 
inertial effects in a perturbative manner into the framework of the classical Saffman- 
Taylor treatment of the viscous fingering instability. 

5.1. Perturbation of the Darcy's law 

Modifications of Darcy's law have already been introduced in S 12.21 We will now consider 
the Euler-Darcy equation and thus a Darcy equation corrected by inertia, in the frame 
of the moving finger. 

One starts from (|2.4|l for the two-dimensional velocity field u{x, y, t) in the laboratory 
frame. In the frame of the moving finger the problem is by definition stationary. Con- 
sidering u{x,y,t) = u'{x — Ut,y) -\- UBx, we obtain for the velocity u'{x,y) in frame of 
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reference of the moving finger: 

(""^S + + Pu'Vu'^ = -Vp - ^u'. (5.1) 

Using the same scaling as before and omitting the * (dimensionless symbols) and the ' 
(finger frame) for the variables, we find: 

Re* ( (/3 - ") ^ + /3m Vm j = -Vp - u. (5.2) 

We assmne that the flow remains a potential flow, i.e. u = V0. This restriction to 
potential flow without vorticity is possible as long as the boundary conditions that apply 
to the finger and the walls are not modified. If this is the case, one can write: 

(f>= -[p+ Re* [{/3 - a)u^ + i/Jw^] +x + est] and A0 = 0. (5.3) 

Considering that along the finger and away from its tip, there are no inertial effects 
(the fluid is at rest in the laboratory frame, and lim^^_oo u = in the finger frame) 
one should choose the constant equal to: Re*{f3/2 — a). 

The mechanical equilibrium of the interface requires the balance of the normal stress 
from both sides, given by ^^-.p^ -j/R (i? > 0), where 7 = {b/W)'^j/{12T]U) = 1/5"^ 
is the dimensionless surface tension and R is the dimensionless radius of curvature. Using 
this, one obtains the following boundary condition for (j) at the interface: 

0r = - [-7/^ + Re* [{f3 -a)u^ + ^(3u^ + (3/2 - a] + x] . (5.4) 

As the normal velocity at the interface it„ is zero in the frame of the moving finger, the 
only remaining velocity component is the tangential one Ut and we can use the notation 
of iMcLean fc Saf fman ( 1981i): u = utet — —q{cos6ex + sinfle^) where q varies from 
(at the tip of the finger) to 1 (at its side) when 9 varies from —n/2 to 0. We can thus 
replace by and —Ux by qcos9. 

As q is mainly given by cos 9, one can write: 

0r = -2; + 7 [1/i? + We* (a - (3/2) sin^ 9] . (5.5) 

N ote that the last term is the Bernoulli correction. In all previous mention ed references 
fsee lGondret fc R,abaudl[l997tlR,uver-(5^l200lHPlouraboue fc Hinchll2002l) {a -(3/2) is 

positive, so this correction has the same sign as the curvature. It also vanishes at the 
sides of the finger. This shows that the effect of the inertial term is very similar to that 
of the capillary forces: the inertial forces should tend to increase the finger width, as was 
indeed observed experimentally. 

Finally, rewriting H5.5|l . it follows that: 

4>r = -X + 7(61) /R, with 7(61) = 7 [1 + We* (a - (3/2)R{9) sin^ 9] , (5.6) 

a form identical to H4.5|l obtained above by considering the modified control parameter 
\/B' with an effective surface tension. 

5.2. Numerical simulations and comparison to experimental data 

Of course the selection of the relative width of the finger can only be found by a sophisti- 
cated singular perturbation analysis. However, the relative finger width can be obtained 
numerically by a modification of the McLean & Saffman (1981) method. We choose this 
method for acomparison with the experimental results and introduce the correction of 
lR,uver-Qiii]l l(200ll) in the numerics using l|5.4|l . We also introduce the effect of the wetting 
film by modifying in the value of the surface tension as done in 1)4. 2|l . 



12 



C. Chevalier, M. Ben Amar, D. Bonn and A. Lindner 




Figure 8. Finger width A as a function of the classical control parameter 1/ B for the geometry 
with &=0.75 mm and W=4 cm and different fluids: numerical simulations (a) and comparison 
between simulations and experimental results, experimental data: », V02; x, V05; □, VfO; 
+ V20; A, VfOO; Numerical simulations: lines. 

If this is done, the numerical simulations confirm the simple argument pointed out 
above and show an increase of the relative finger width compared to the classical results 
(see figure IHt) in agreement with the experiments. The observations from the numerical 
results are: 

• The inertial effects are stronger (i.e. appear for a smaller critical 1/ B) for less viscous 
fluids as well as for larger cell thickness or for smaller cell width. 

• The minima of the relative width as a function of We* are around a unique value 
of M^e*, however the numerical value of We*^ differs between the experiments (^15) and 
the simulations («2). 

When comparing the simulations and the experiments (see figure|SlD) it turns out that 
they are in qualitative agreement but that the results are not identical. 

When inertial effects are present we can characterize these effects by estimating a 
critical value of the control parameter 1/ Be and a critical relative width Ac at the 
minimum. It turns out that the numerics provide a rather good estimate for Ac but fall 
to give a correct value for 1/ Be which is found to be smaller in the numerical simulations 
than in the experiments. Another significant difference is that the increase of the X-1/ B 
curve is observed to be stronger in the simulations. 

However, even for the case where inertia can be neglected (for most viscous fiuid), there 
is a small but significant discrepancy between the numerical simulations and the experi- 
mental data. We believe this is due to the fact that our way of correcting for the wetting 
film is too simple. In consequence, it is clear that we can not expect perfect agreement 
between experimental data and simulations when adding corrections due to inertia. To 
quantitatively account for the experimental results, one has certainly to go back to the 
three-dimensional effects of the experiment which are expected to be important close 
to the finger tip over a length scale of order b. This is due to the existence of a three- 
dimensional structure of the fiow which ca nnot be ignored in the vici nity of the finger: a 
film ex ists between the plate and the finger. IPark fc Homsvl l)l984(l and lReinelt fc SaffmarJ 
l|l985|) have shown that it is nevertheless possible to reduce the prob lem to two dimen- 
sions when modifying the boundary conditions on the finger (see also lBen Amar fc Ricd 
1200211 . However, they have also shown that reduction to two-dimensional is only possi- 
ble if the parameter We = pU^b/^ is small. For our problem, this is not the case and 
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therefore a complete set of the corrected boundary conditions must be deduced from the 
full three-dimensional theory of Park, Homsy and Reinelt incorporating inertial effects 
in order to resolve the problem, which is largely beyond the scope of this paper. 

6. Summary and conclusion 

Wc have investigated the effect of inertia on the Saffman-Taylor instability. Inertial 
effects are found to become important for high Reynolds numbers and thus for fluids 
of low viscosity and for large plate spacing of the Hele-Shaw cell. For these situations 
one observes, upon increasing the velocity, first a classical regime with a decrease of 
the relative finger width and then a second and new regime in which the finger width 
increases. This second regime is due to the importance of inertia. 

We introduced a modified Weber number We* which allowed us to explain the crossover 
between the two regimes. The transition is thus given by a critical modified Weber number 
We*. Below We*, the classical regime of decreasing finger width is of course governed 
by the classical control parameter 1/S, which is a modified capillary number. The finger 
width in this regime is thus given by the balance between capillary forces, which tend 
to widen the finger, and viscous forces, which tend to narrow the finger. With increasing 
velocity the viscous forces dominate over the capillary forces and one observes a narrowing 
of the finger. For the second regime, above We*, one observes on the contrary an increase 
of the finger width with increasing velocity. In this case the finger width is governed by 
a modified Reynolds number Re* and thus by the balance between viscous forces and 
inertia. It turnsout that inertial forces tend to widen the finger. With increasing velocity 
inertia dominates the viscous forces and one consequently observes a widening of the 
fingers. 

We have also shown that we can define a new control parameter l/B\ which takes the 
corrections due to inertia into account. This parameter tends towards 1/B for low We* 
and is proportional to 1/ Re* for large We* . When plotting our data as a function of this 
empirical parameter they collapse onto a single universal curve which corresponds to the 
results for the finger width obtained for the classical Saffman-Taylor instability. 

By only taking into account a modification of Darcy's law, some simple arguments and 
numerical simulations confirm all of these observations. However, the agreement between 
numerics and experiments is only qualitative. We believe this is due to the fact that the 
problem is certainly three-dimensional and one must consider the full three-dimensional 
theory of Park, Homsy and Reinelt incorporating inertial effects. 

We thank Eric Clement, Mike Shelley and Laurent Limat for usefull discussions and 
Jose Lanuza for valuable help with the experimental set-up. 
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